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SUMMARY 

A theoretical analysis of absorbers whose radioactivity is 
time-dependent is presented. Matrix methods are employed as a tool 
in the analytical determination of the intenstiy of radioactivity 
and the amount of heat generated in any portion of a thick absorber. 
The method is applied to the following three cases: (l) plane 

source of monochromatic radiation of a single type at normal inci- 
dence to a plane absorber, (2) plane source of polychromatic radia- 
tion of a single type at normal incidence to a plane absorber, and 
(3) plane source of polychromatic radiation of several types at 
normal incidence to a plane absorber. Illustrative examples are 
included. 

Experimental data on -the absorption of gamma rays and neutrons 
and on the rates of conversion of these radiations to thermal energy 
are required for utilization of this method. 

INTRODUCTION 

A mathematical analysis of the interaction of nuclear radia- 
tions with matter requires simultaneous consideration of a large 
number of factors. Fortunately, many of these factors obey the 
same mathematical laws. Hie compact notation inherent in matrix 
methods makes them ideal for problems of this nature. 

The matrix method was first used to calculate the transmission 
and reflection of thermal radiations from the surfaces of a number 
of thin foils (reference 1). A similar method has been applied to 
nuclear radiations and developed for solutions of cases of poly- 
chromatic incident radiation, neutron and gamma-ray combinations, 
radiation not normal to the face of the absorber, and time- 
independent radioactivity in the shield itself (reference 2). 
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Inasmuch as a great many artificially produced radioactive 
substances have very short half-lives , that is, the radioactivity 
of the substances decays rapidly, the extension of the method to 
include time-dependent cases was considered useful in widening its 
applicability. Thus, an extension of the method was developed at 
the NACA Lewis laboratory to include the following time-dependent 
cases: 

(1) The incident radiation is monochromatic. The radiation is 
normal to the absorber. Part of the radiation absorbed at each 
station of the plane absorber is transformed to heat and part to 
radioactivity in the absorber. The absorber emits radiation of the 
same type and wavelength as the incident radiation. The radio- 
activity produced in each station of the absorber is entirely 
emitted. 

(2) An extension of case (1) to include polychromatic radiation 

(3) An extension of case (2) to include mixtures of gamma and 
neutron radiation 

As used herein, the expression "type of radiation" refers to 
the fundamental particles emitted by radioactive elements. Gamma 
radiation is one type of radiation; neutron radiation is another 
type. 

Case (1), which assumes no degradation of energy, is presented 
for the sole purpose of illustrating the use of the method with a 
minimum of cumbersome notations and computations. 

Experimentally verifiable data on all the constants used must 
be supplied for proper utilization of the method. Such constants 
include absorption coefficients, reflection coefficients, trans- 
mission coefficients, rates of conversion of energy of nuclear 
radiations to thermal energy, and radioactive decay constants. For 
gamma rays, most of the required absorption data are already avail- 
able or can be computed; but for neutrons, the amount of available 
data, although large, is still much smaller than is to be desired. 


ANALYSIS 

This analysis may be conveniently* divided into three separate 
cases. Each of these cases, together with illustrative examples, 
are subsequently described. All symbols are defined in appendix A. 
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Case I 

The assumptions used in case I are: 

(1) Stations of absorber are parallel plane surfaces. 

(2) Incident radiation is monochromatic. 

(3) Radiation is normal to absorbers. 

(4) No degradation of energy occurs. 

(5) Part of energy of radiation absorbed is transformed to 
energy of radioactivity and part to thermal energy. 

(6) Absorber emits radiation of same type and energy as 
incident radiation. 

(7) Radioactivity produced in each station of absorber is 
entirely emitted. 

Method. - Diagram (a) is included to clarify the discussion 
in the following paragraph. 



The absorber may be 
considered as consisting 
of n stations normal to 
the path of the radiation 

The i^* 1 station of the 
absorber is denoted by 
Sj^. Radiation Ij^ from 
the preceding station 
S 1 _-^ Is incident upon 

one surface of 
while radiation R i+1 
from the subsequent sta- 
tion S i+ ^ impinges 
upon the opposite sur- 
face of From S ^ 

emerges radiation 


which is incident upon Sj_ + ]_, and radiation R^, which is incident 


upon 
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In the treatment presented herein, the following quantities 
are assumed to he known time- independent constants: 


= X £ 
^ S i^T = Tf 


x i 
**i 

*i 

^i 


radiation power incident upon absorber from source 

radiation power incident upon absorber from surround- 
ings outside absorber 

one-half of initial power of radioactivity in each 
station of absorber 

initial thermal power being emitted from each sta- 
tion of ' absorber 

power-transmission coefficient 

power back-scattering coefficient 

rate of conversion of energy of radiation absorbed 
to energy of radioactivity 

radioactive decay constant 

rate of conversion of energy of nuclear radiation 
absorbed to thermal energy for each station of 
absorber 


The problem is to find 1^ and Rj[ in terms of known quan- 
tities. The following procedure is used to solve the problem: 

All Ij_ and R^ are expressed in terms of I^_, R n+ -^, tj, 
r^, and Pj_. These expressions are then substituted in the dif- 
ferential equations for Pj_ and the resulting set of n simulta- 
neous linear differential equations of the first order are solved 
for P 1# This procedure will yield P-j_ in terms of known quan- 
tities; when these values of Pj^ are substituted in the expres- 
sions for Ii and R i; the original problem will be solved. 

If these expressions for Ij_ and R^ are substituted into the 
equation expressing the conversion of radioactive energy to thermal 
energy, Hj[ can be found. 

For the i"** station of the absorber Sj_ (diagram (a)), the 
following relations hold: 
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I i+1 = t i I i + r i R i+1 + p i (!) 

R i = r i I i + t i R i+1 + p i (2) 

% = hi (Ii + E i+1 ) (3) 

dP-i 

— = n ± (Ii + R i+1 ) - ^Pi (4) 


The derivation of equation (4) is given in appendix B. 

Equation (4) is the basic feature of the treatment of time- 
dependency in this paper. The absence of additional terms in this 
equation (and also in equation (3)) is due to the simplifying 
assumption made that the radioactivity produced in each station of 
the absorber is entirely emitted. In practice this assumption is 
effectively true in most shielding problems. 

The use of equations (1) and (2) in expressing Ij_ and R^ 
in terms of and known quantities is explained: 

Equations (1) and (2) may be rewritten in the following 
manner: 



In matrix notation, equations (5) and (6) may be expressed as 



If the square matrix of constants is denoted by [t^ ] and the 
column matrix of constants by equation (7) may be simplified 

to 
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x i 

R i 


- M 


± i+l 

R i+1 


+ ^i|. p i 


( 8 ) 


Beginning with the expression for 


tion of expressions for subsequent 


R i 


, continued substitu- 


matrices will result in 


I 1 

-(tt [4 

P n+1 

_ R 1_ 

W / 

_ R n+ 1 _ 


+ s(H n) n ^ 


O) 


where jjr 0 [| is defined as the identity matrix [e]. 
Similarly, equation (7) can be written as 


li+l 


R i+1 


r 2 ^ r 
t 4 - r-i 


'i t 


t i 


E, 


l-SL 

*i 


*i 


p i 


(7a) 


The square matrix of constants 
the column matrix of constants is 
tion (7a) may be obtained directly 
algebra . 


If the analogy is extended, 


I n+1 

= fnW- 1 ) 

V 

R n+1 

11 

P 

1 

1 

X 

R 1 


is actually [Ti] - "*", whereas 
[Ti]" [Qi] so that equa- 
from equation (7) by matrix 
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Although in the following discussion equations of the form of 
(9a) will he arbitrarily employed to obtain the t Iq and Rq in the 
desired forms, the use of equation (9) leads to the same results. 

If matrix notation is temporarily abandoned, equation (9a) is 
written as two separate equations of the form 

n 

I n+1 = a l z l + a 2®l + ^ Vj (10) 

\tl = a 'l I l + a ’2»l + (ID 

j— l 

Inasmuch as Iq and R n+ q are constants, equations (10) and. 
(11) may be seen to be in the forms desired. Similar expressions 
for other Iq and Rq may be found as follows: 

The following equation can be written from equation (7a): 

12 1 [ Ix l 

* h]- 1 - [Ti]-1 [Qi] Pi (12) 

r 2 r i 

Writing out the two equations represented by equation (12) and 
substituting equation (ll) for Rq results in expressions for I2 
and Rg in terms of known quantities and Pq. These expressions 
may then be used to find I3 and R3, which in turn enables the 
determination of I4 and R4 and so forth until all Iq's and Rq's 
are in the desired forms. 

After the expressions for Iq and Rq, which are linear 
combinations of Pq have been obtained, the next step is to sub- 
stitute these expressions into forms of equation (4). The result- 
ing set of n differential equations is of the form 
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~ = C 1 I 1 + d l R n+l + ^ e l,j P 0 


cLPo n 

£ . T \ 


— = cglx + <i 2 R n+1 + Z) ®2, j p j 

J — J- 


a ii 

It" 15 c i I l + *l*n+l + ^ e i, j p j 

j — 1 


df = c n p l + ^1 + Z, ®n, j p j 

/ 

where all c's, d.'s, and e's are arithmetic combinations of t i# 
ri, Xi, and (ii, so that all quantities on the right-hand side 
of equations (13) are constants except the P's. 

Equations (13) may be written as the single matrix equation 

^ - M +[aJM U< 


where [ P] is the n- rowed column matrix f P^ 
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[A] is an n-rowed square matrix, all of whose elements are con- 
stants, and [ B] is an n-rowed column matrix all of whose elements 
are linear combinations of I]_ and R n+1 and therefore constants. 

The solution to equation (14) is obtained in the same manner 
- as if all its terms were scailars. 

Thus, if P = P° at t = t q , as long as [a] is non- 
singular, the solution to equation (14) is 

[P)-= [A]’ 1 (eW(T-To) ([ A ][p0] + [B]) - IB]} (15) 

(See appendix C for explanation of e^-K T " T o) . ) 
where C P^] is the n-rowed column matrix 


After P^'s have been found in terms of known quantities, 
these P^s may then be substituted into expressions of /the form 
of equations (10) and (ll) to find 1.^ or Rj_. If these expres- 
sions are substituted into equation (3), may be found. 

Example. - Assume that 

(1) The absorber is three stations in thickness. 

(2) The incident radiation is monochromatic. 

(3) The radiation emitted by the absorber is of the same type 
and energy as the incident radiation. 



\ 
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(4) I 1 = 0.1 roentgen per hour 
. r 4 = 0.01 roentgen per hour 
t^ — 0.8 


0.5 


t 3 = 0.8 
r^ = 0.1 
r 2 = 0.4 

r 3 = 0-1 

hi = hg = h 3 = 5 x 10 


,-2 


V -1 = ^2 = ^3 


= 5 x 10" 7 seconds -1 


as lO" 5 seconds -1 
Ag = 2 X 10“ 5 seconds” 1 


X 3 = 5 X 10” 5 seconds 


-1 


(5) At t = 0, P° = 


and H x = Hg = H 3 = 0.00. 


The problem is to determine I4, and H2 at t = 10 

seconds . 

Solution. - If the given data are substituted into equa- 
tion (12) there results 


,5 - 
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Simplifying equation (16) and solving for I* and R^ will result 


in 

I 4 = 0.0387 + 0.436 P x + 0.880 P g + 1.356 P 3 (17) 

E 1 = 0.0420 + 1.356 P x + 0.880 P 2 + 0.435 P 3 (18) 

Similarly, 

1 2 = 0.0840 + 1.0445 P x + 0.110 P 2 + 0.0544 P 3 (19) 

R 2 = 0.040 + 0.445 P x + 1.10 P 2 + 0.544 P 3 (20) 

1 3 = 0.047 + 0.544 P L + 1.099 P 2 + 0.445 P 3 (21) 

E 3 = 0.0127 + 0.0544 P x + 0.1099 P 2 + 1.0445 P 3 (22) 


From equation (4), the rates of change of P-^ P 2 , and P 3 are 

given hy 
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dP-, 7 c 1 

_± = 5 X 10“' (0.10 + Eg) " 10" 5 P]_ 

dPo 

= 5 X 10 (ig + E 3 ) - 2 X 10” 5 Pg > 
dP 3 _ 

^ =5X10'' (I 3 + 0.01) - 5 X 10“ 5 P 3 


(23) 


where Ig, Eg, I 3 , and E 3 are given by equations (19), (20), 
(21), and (22), respectively. 

Substituting these expressions into equation (23) and con- 
verting to matrix notation results in 


d[p] 

dT 


= 10 


-6 


-9.778 

0.549 

0.272 


0.070 

.549 

-19.890 

.549 

[p] + 10” 6 

.048 , 

.272 

.549 

-49.778 


.0285 


(24) 


Applying equation (22) and writing out the equations for each 
P^ results in 

P x * - 7.42X10" 3 e^ lT +1.22X10” 4 + 3.48X10” 6 e^ T + 7 .29xl0” 3 

Pg = - 4. 06x10” 4 e^ lT - -2 . 24x10” 3 e^ T + 9.87X10” 6 e^ T + 2.63xl0” 3 ^ 

P 3 - - 5. 62x10” 5 e^ T - 4.00X10” 5 e^ T - 5.46X10” 4 e^ T + 6.42X10” 4 


where 

= - 9.746 x 10” 6 

1 2 = - 19.910 X 10" 6 

1 3 = - 49.790 X 10” 6 


( 25 ) 
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Substituting equation (25) into equations (17) an d (18) will 
give I 4 and E^. Thus at T = 10 5 seconds, 

I 4 = 0.0434 roentgen per hour (26) 

E^ = 0.0527 roentgen per hour (27) 

By the use of equation (25), first I 2 and B3 and then 1 2 

may "be found. If the given data and the newly found expressions 
for I 2 and B3 are substituted into equation (3), 

H 2 == h 2 (I 2 + E 3 ) 

= 5 X 10" 2 (0.0850 + 0.0138) 

= 4.94 x 10" 4 roentgen per hour 

Case II 

The assumptions used in case II are: 

(1) Stations of absorber are parallel plane surfaces. 

(2) Eadiation is polychromatic and normal to absorber. 

(3) Energy i3 degraded. 

(4) Part of energy absorbed is transformed to energy of radio- 
activity and part to thermal energy. 

(5) Absorber emits radiation of same type as incident radiation. 

(6) Eadioactivity produced in each station of absorber is 
entirely emitted. 

Method. - In order to make the method of case I amenable to an 
analysis of polychromatic radiation, the terms are defined as 
follows: 

The energy spectrum of the radiation, which includes the 
spectrum of the radioactivity in the absorber, is divided into a 
finite number m of energy bands. Each band is assigned a single 
energy, which represents all energies within the band; for example, 
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if the limits of a hand are 0.50 million electron volt (Mev) and 
1.49 Mev, and if 1.00 Mev represents all energies in the hand, then 
energies such as 0,51, 0.70, and 1.20'Mev would all he treated as 
energies of exactly 1.00 Mev. By varying the number of hands, any 
desired degree of accuracy may he attained. 

n 

From the preceding paragraph, Ij[ = vhere ij^j is 

^“1 

the power of radiation within the energy hand incident upon 

the i"*^ station of the ahsorher. Each 
in the form of an m-rowed column matrix 


I ± may also he expressed 


Ci ±3 = 


L i,i 




Similarly, Ri and Pj_ may he expressed as the m-rowed 


column matrices 


_p i,r 

• 

and 

"Pi,r 

« 

* 


• 

• 

• 

p i,m 


• 

• 

p i,m 


, respectively. 


With these considerations, relations analogous to equations (1) 
through (4) may he defined as follows: 


Pi+l] = C'fciD C-t±] + M [R i+1 ] + [Pi] 

[Bj] = [rj][li] + [tj] [Ri + ]_]’ + [Pj] 


(28) 


(29) 
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= CnilCli + E i+1 ] - [X.J] [P^ (30) 


[%] = [hi] [ii + E 1+1 ] (31) 

where [tj.], [ri] , [pi],[^j], and [hj] are' m** 1 order square 
matrices, whose elements tj,k, rj k , k> and bj i> 

respectively, denote (l) the fraction of incident radiation in 
energy hand k transmitted to energy band J; (2) the fraction of 
incident radiation in energy band k back-reflected to energy 
band j; (3) the rate of absorption of radiation power from energy 
band k to become power of radioactivity in energy band j; 

(4) the decay constants of radioactivity in the k tb energy band; 
and (5) the rate at which nuclear radiation in the energy band j 
is transfoimed to thermal energy. 


Before continuing the analysis, the forms of the aforementioned 
matrices may be simplified. For this purpose, the convention will 
be adopted of dividing the energy spectrum into bands of increasing 
radiation energy so that the band of lowest radiation energy will 
be the first band while the m^b band will be the highest radia- 
tion energy in the spectrum. By so doing, [t A ] and [r*] can be 

reduced to triangular matrices. Inasmuch as transfer of radiation 
from any energy band to a higher one can occur only after a nuclear 
reaction, and because the results of all nuclear reactions will 
appear as p^^ or H^, there will be no elements of [ t 1 ] or 
[r^] denoting an upgrading of energy. 

This observation is general and applies to any type of nuclear 
radiation. However, further reduction of [tj^ ] or [ r^] to 
diagonal matrices is prevented if allowances are to be made for 
loss of energy due to either elastic or Compton back- scattering. 

No great simplification can be made for [ p^] . However, at 

present, there appears to be no experimental evidence of any radia- 
tion that excites a radioactivity of both the same type and energy 
as itself. This fact would serve to put zeros in the diagonal ele- 
ments of [pj_]. 

By definition, [A ^] and [ h^] are diagonal matrices. 

In order to continue the analysis, equations (28) and (29) are 
solved for 1^ and in terms of I^ + ^: 
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[Ij = [td' 1 [I 1+ i] - [td* 1 [rdOBl+d - [td* 1 [Pd (32) 
[Ed * bdttd' 1 [I i+ d + ([td - [riDCtd-Vi])^] ' 

+ ([E] - [ r d [ t d ~ 1 J[ p l] (33) 

where [e] is defined as the identity matrix of order m. 

Inasmuch as all terms in equations (32) and (33) are matrices, 
the order of multiplication must he carefully preserved. The 
assumptions made here imply the existence of the inverses of var- 
ious matrices. Whenever a solution to the physical problem is 
possible, all matrix inverses as used in this paper will exist. 

In matrix notation, equations (32) and (33) may be written 


li 


^d* 1 

- [‘d* 1 [ r d 


I i+1 

A. 


i 

i — i 

►t-i 

w 

H* 

1 — I 

1 

[tj] - [rd [td* 1 [rd 


E i+1 


- IHT 1 

[«]- [rd&d' 1 


(34) 


If the notation is simplified, 


V 


I i+1 


-[Tj 


Bi 


E i+1 


[QjJCPi] 


(35) 


where [ ] is a 2m-rowed square matrix and 

rectangular matrix. 

*1 

Beginning with the expression for 



tution of subsequent I i and B^ yields 


[Q^] is a 2m-by-m 
, continued substi- 
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V 


I n+l 

_ K 1_ 

/ 

_ E n+l 


+ 




J=1 



fy3l>j2 


Both J and k must he taken in the order indicated. 


(36) 


There are 2m equations in 2m unknowns so that the i n+ i 
and the p, may he found in teiros of the p. . and known quan- 

J- 1 j J 

tities. 


From this point on, the method is the same as the procedure 
in case I beginning with equation (12). 

Example. - Assume that: 

(1) The ahsorher is three stations in thickness. 

(2) The incident radiation is composed of two monochromatic 
radiations of equal intensity. 


(3) There is no hack- scattering to the absorber. 

(4) The radiation emitted by the absorber is of the same type, 
but of a different energy, than the incident radiation. 


(5) 




1 

1 

0 


R 4 


(P° 


i 



0 

0 

0 


(Hi)T 0 =0 « 


0 

0 

0 


CtJ 

Ct 2 ] 


0.50 

0.00 

0.00 

.01 

.60 

.00 

.01 

.01 

.50 

0.60 

0.00 

0.00 

.01 

.50 

.00 

.01 

.01 

.50 
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[t 3 ] = 




[r 2 ] = 


l>*] = 


0.50 

0.00 0 

.00 

.01 

.50 

.00 

_ .01 

.01 

.60_ 

'0.40 

0.00 0 

.00" 

.01 

.30 

.00 

_ .01 

.01, 

.30 

'0.30 

0.00 0 

.00" 

.01 

.40 

.00 

.01 

.01 

.40_ 

"0.40 

0.00 c 

1.00 

.01 

.40 

.00 

. * 01 

..01 

,30_ 


"0 0 0" 


10“ 2 

0 0 0 
111 

sc 


,-1 


[**2 3-10 


,-2 


0 

0 

0 


second 


-1 


v*2 


0 

0 

1 


[1*3 = 10 

C^3 - s C^3 - 10-1 

Ol3 - [^2^ = [^3^ = 10 2 


second 


-1 


0 

0 

1 

0 

0 

1 


0 

0 

1 

0 

1 

0 


0 

0 

1 

1 

0 

0 


second' 


-1 


The problem is to determine I4 and at t = 1000 seconds 

If the above data are substituted into equation (36), 


= [TJ[T 2][T 3 ] 


*4 

R 4 


CQiDCpJ + Ct 1 :[q 2 3 Cp 2 D 


[TjCTgHQslCPj] 


(37) 
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where 


M- 


[Qj- 


[T 2 ] = 


2.000 

0.000 

0.000 -C 

o 

o 

00 

• 

0.000 0.000 

-.033 

1.67 

.000 - 

.00335 -.500 .000 

-.039 

-.033 

2.000 - 

.00395 -.010 -.600 

.800 

.000 

.000 

.180 

.000 .000 

.01 

.500 

.000 

.001 

.450 .000 

.00788 

.00668 

.600 

.00078 .002 , 320 j 

- 2.000 

0.000 

o . ooo ' 




.033 

-1.67 

.000 




.039 

.033 

- 2.000 




.200 

.000 

.000 




-.010 

.500 

.000 



■'* 

-.00788 

-.00668 .400 




- 1.67 

0.000 

0.000 -0 

500 

0.000 o . ooo " 


-.033 

2.000 

.000 - 

010 

-.800 .000 


-.033 

.040 

2.000 - 

010 

-.036 -.800 


.500 

.000 

.000 

450 

.000 .000 


.00335 

.800 

.000 

obi 

.180 .000 


.00326 

.004 

.800 

001 

-.012 .180 
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[ t 3 3 = 


C%3- 


1.67 

0.000 

o . ooo " 



.033 

- 2.000 

.000 



.033 

-.040 

- 2.000 



.500 

.000 

.000 



-.00335 

.200 

.000 



-.00326 

-.004 

.200 



2.000 

0.000 

0.000 

- 0.800 

, 0.000 

-.040 

2.000 

.000 

-.004 

-.800 

-.0333 

- 0.0333 

1.667 

-.00326 

-.00335 

.800 

.000 • 

.000 

.180 

.000 

.004 

.800 

.000 

.0004 

.180 

.00979 

.010 

.500 

.001 

.001 


- 2.000 

0.000 

0.000 

.040 

- 2.000 

.000 

.033 

.033 

- 1.67 

.200 

.000 

.000 

-.004 

.200 

.000 

-.00979 

-.010 

.500 


If the indicated operations are performed, 


0.000 

.000 

-.500 

.000 

.000 

.450 
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I 

• 

CD 

0 

0 

-2.60 

0 

0 




.353 

4.73 

0 

.079 

-2.60 

0 


V 


-.390 

-.430 

4.668 

.0929 

.114 

-2.60 


V 

E l_ 


2.590 

0 

0 

-1.20 

0 

0 


*4_ 



-.0589 

2.46 

0 

-.009 

-1.15 

0 





-.0755 

.0806 

2.23 

-.0022 

.0002 

-1.03 



+ 


-2 

0 

-i 

0 


-3.73 

0 

0 

.033 

-1.67 

0 


.111 

-3.43 

0 

.039 

.033 

-2 

[Pi] + 

.128 

.145 

-4.00 

.2 

0 

0 

-1.24 

0 

0 

1 

• 

o 

H 

.5 

0 


-.001 

-.91 

0 

-.00788 

-.00688 

.4 


.0071 

-.011 

-1.24 


[P 2 3 


-6.14 

.351 

.388 

-2.91 

.049 

.067 


0 

-6.15 

.449 

0 

-2.78 


0 

0 


-6.67 

0 



0 


.073 -2.83 


(38) 


Substituting the values for I-|_ and R 4 into equation (38) 
and solving for I 4 and R]_ yields: 



— 0.209 + 0.418p^^ + 0.781p 2 -^ + 1.285p2^ 
i 42 =0.227 + 0.0242p 11 + 0.353p 12 + 0.0348p 21 + 0.726p 22 + 0.0217p 31 
+ 1.301p_„ 
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09366 -.09352 -.09288 .01642 .01503 .,01430 .00906 .00735 .00712 
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The solution to equation (15) is thus: 

P 11 = P 12 = P 21 = P 22 = P 31 = P 32 = ° 

p = -0. 1115e -0 * 1000 T - 0.03639e“ 0,1090T ' 

r 13 

- 0.1779e"°* 068:L8T + 0.3257 

p 23 = 0 . 0000468e”° * 1000T + 0.05914e" 0, 10901 ► ^ 

- 0.2187e" 0,06818T + 0.1595 

P 33 = 0 • 1113e "° * 1000t - 0.03626e"°* 1090T 

- 0.1779e”°*° 6818T + 0.1030 j 

If these expressions are substituted into equation (39), at 
t = 1000 seconds 


141 = 0.209 roentgen per hour 
*42 = ° .227 roentgen per hour 
j"43 = 0.329 roentgen per hour 


After Eg has ,been found in a similar manner, the resulting 
expressions may he substituted into equation (31) from which 
is found to be H-j_ = 0.0326 roentgen per hour. 


Case III 

The assumptions used in case III are: 

(1) Eadiation consists of polychromatic neutron and gamma 
radiation. 

t 

(2) Eadiat'ion is normal to absorber. 

(3) Part of energy absorbed is transformed to energy of radio- 
activity and part to thermal energy. 
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(4) Absorber emits radiations of same types as incident 
radiation. 

, (5) Radioactivity produced in each station of absorber is 
entirely emitted. 

Method. - The method differs little from that suggested in 
case II. Each type of radiation is dealt with separately, whenever 
possible. The energy spectrum of the gamma radiation is divided 
into m energy bands and that of the neutron radiation into q 
energy bands. These energy bands have the same denotations as in 
case II. In the analysis that follows, symbols for gamma -radiation 
power will be signified by the superscript y and those for 
neutron-radiation power by the superscript n. As may be readily 
seen from the following description, the method is immediately 
applicable to any number of 'different types of radiation. 


Si 



The following relations hold for this case: 

[I r 1+ l] - [t 7 i]tl 7 i] + O 7 !] [K r 1+ l] +[P 7 l] (41) 

[K 7 !] - [r 7 i][I 7 i] + Ct 7 1 ][E 7 1+1 ] + [P 7 i 3 (42) 


= 0 7 i3 1 i 7 i 4 R 7 i + 1 3 + O n ’ 7 i3 [I n i + R n i+ i3 - C ^ 7 i3C? 7 i3 


( 43 ) 
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P n i+x3 = P n J[i n i3 +Cr“iJ_- <«) 

[R“l] = [r^C I“t ] ♦ [t n i]CB n i + l] + [P"t] < 45 > 

iHji = Cv n i][i n i + B"i + i] + O r '“de r i + K 7 i +1 ] -C^ n i ][P n i] 

4t 

(46) 

[%] = [h^][l7. + H7 1+1 ] +[ h n i][I n i + R n i+1 ] ( 4? ) 


where [t 7 ^], [r^], [n 7 *], [X r ± ] , and [h r ± ] have the same mean- 
ings and are in fact identical with [tj], [rj], [jiJ, [X*], and 
[ hj , respectively, in case II; the corresponding constants for 
neutron radiation have analogous meanings: [ |a n > 7 ' i ] is an m-by-q 

rectangular matrix whose elements p n . 7 ^ denote the rate at which 

neutron radiation of energy j is converted to radioactive power 

of gamma radiation of energy k, and correspondingly, [ ] is 

a q-by-m rectangular matrix whose elements u denote the rate 

'u“v 

at which gamma radiation of energy u is converted to radioactive 
power of neutron radiation of energy v. 

The set of equations (41) and (42) is solved to obtain all 
[I 7, ] and [R 7 ^] in terms of known quantities and [ P?] by the 
same procedure as of case II. Then equations (44) and (45) are 
treated similarly so that ,[I n .] and [ R n ^] are obtained as expres- 
sions consisting of linear combinations of known quantities and [P 1 ^] 

✓ 

The resulting expressions are substituted into the set of equa- 
tions (43) and (46), which is combined to form the single matrix 
equation 

= [A] CP] + [B] 

where now [?] is the n(m+q) -rowed column matrix 
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/V 

and [A] and [B] are n(m+q) -rowed square and column matrices 
of constants, respectively. 

After this point is reached, the methods of case I may he 
applied to complete the solution. Although physical problems might 
occur where the method of case I will not yield unique results, 
these cases are very rare. 

The preceding method is immediately applicable to an analysis 
of radiation consisting of any number of different types of radia- 
tion for which the required experimental data are available. 


DISCUSSION 

Advantages of the Method 

The advantages of the foregoing method are twofold. First, 
the inclusion of time-dependent factors greatly increases the num- 
ber of problems to which the method may be applied. Almost all 
chemical elements have some radioactive isotopes. Inasmuch as many 
radioactive substances have short half-lives, this method is a 
closer approximation to the physical problem than a method that 
considers the radioactivity to be a time- independent constant. 

Another advantage is that all the P-^'s may be obtained with 
approximately the same amount of calculations required to obtain a 
single without the use of matrices. If matrix methods are 

not used, the usual procedure is to obtain an n ti] order equation 
in one Pj^ by substitution, solve it by differential operators, 
and then substitute this solution into the expressions for each of 
the other P-^'s in terms of the one just obtained. 
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If computing machines are available, the calculations of [A]”'*' 

and e CA]T do not take long and the matrix method may be readily 
employed to the fullest advantage. 

An example of a possible application of this method is: 

Assume that two prospective absorbing materials are available 
for diminishing the intensity of a certain radiation. One material 
is an excellent absorber but rapidly becomes radioactive when sub- 
jected to the given radiation. The other material, although 
inferior to the first in its ability to absorb the radiation, lacks 
the undesirable property of becoming radioactive. Which absorber 
should be chosen? 

Obviously, there are situations in which the second material 
would be preferable and other situations in which the first mate- 
rial would be preferable. The use of this method could shorten the 
time spent in making choices of this nature. 


Limitations of Method 

The main limitation of the method is the amount of available 
experimental data on the absorption of gamma rays and neutrons and 
on the rates of conversion of these radiations to secondary radia- 
tions and to thermal energy. 


CONCLUSION 

A matrix method for determining the effectiveness of radio- 
active absorbers has been set up; theoretically the method is 
immediately applicable for solution of absorber problems. Ade- 
quate experimental data on the absorption of radiation and rates 
of conversion to secondary radiation and thermal energy are, how- 
ever, required for utilization of this method. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, February 4, 1949. 
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APPENDIX A 


SYMBOLS 

The following symbols are used in this report: 


1 — 1 

> 
t — 1 

v» 

1 — 1 
a 

matrix coefficients 

a, a ' ,b,b * ,c,d,e 

scalar coefficients 

l 1 

w 
1 1 

identity matrix 

Hi 

thermal power generated in i^h station of 
absorber 

1 1 

•H 

w 
1 — 1 

matrix of thermal power generated in i^H sta- 
tion of absorber 

h 

rate of conversion of energy of nuclear radiation 
absorbed to thermal energy 

[hi] 

matrix of rate of conversion of energy of nuclear 
radiation absorbed to thermal energy 

*i 

power of radiation from (i-l)^H station 
incident on i™ 1 station 

< 1 

•H 

M 

i — i 

matrix of power of radiation from (i-l) t53 sta- 
tion incident on i^* 1 station 

*1+1 

power of radiation from i til 'station incident on 
(i+l) th station 

Cx i+ i3 

li. 

matrix of power of radiation from i 1 station 
incident on (i+l)^ 1 station . 


power of radiation in energy band incident 

on i^H station from (i+l)^* 1 station 

n 

number of stations in absorber 

[p] 

matrix of power of radiation from radioactivity 


/ 



30 

[P°] 


[%] 

% 


[Ei] 


' 1+1 


[ E i + i3 


[*±] 

s i 

[t 0 3 

[TJ ' 

t 
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matrix of initial power of radiation from radio- 
activity 

) 

one-half of power of radiation from radioactivity 
in i^* 1 station 

matrix element of power of radiation from radio- 
activity in the station 

rectangular matrix coefficient 

power of radiation from i tla station incident on 
(i-l)^* 1 station 

matrix of power of radiation from i^* 1 station 
incident on (i-l) tla station 

+Vi 

power of radiation from (i+1) station inci- 
dent on i tb station 

^ \ 

.LI- 

matrix of power of radiation from (1+1)^ sta- 
tion incident on i” n station 

power hack- scattering coefficient 

matrix of power hack- scattering coefficient 

J_T_ 

i n station of ahsorher 
identity matrix [E] 
square matrix coefficient 
power-transmission coefficient 


[t^] matrix of power- transmission coefficient 

A decay constant for radioactivity, seconds”^- 

[X i ] matrix of decay constant for radioactivity, 

seconds ^ 

' rate of conversion of energy absorbed in to 

energy of radioactivity in 
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TT 

p i,j 


matrix of rate of conversion of energy of radia- 
tion absorbed in Sj_ to energy of radioactivity 
in Si 

symbol denoting that product of terms following it 
are to be taken as indicated by indices 

power of radiation in energy band emerging 

from i^k station 


T 


time 


T 0 

Subscripts: 

T 

T o 

Superscripts : 
n 


initial time 


time 

initial time 


neutron radiation 


7 


gamma radiation 
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APPENDIX B 


DERIVATION OF EQUATION FOR POWER OF RADIOACTIVITY 
IN i th STATION 


A station of an absorber is assumed to consist of a very- 

large number of stable atoms and some radioactive atoms and all 
radioactivity produced in this station is assumed to be emitted. 
The rate of change of the number of radioactive atoms N^ in the 
absorber is given by 


dN i 

dT 


Ws,R 


*i N i 


(Bl) 


where a.^ equals the fraction of the radiation energy incident 
upon the i*“ station that is absorbed in forming radioactive atoms 
and (E^) v equals the energy required to change a stable atom to 

a radioactive one. v , 

Now 


2P i = < E i>E,s < B2 ) 

where C%)rj S is the energy emitted when a radioactive atom 

changes to a stable one. If equation (B2) is differentiated with 
respect to time, 

- dP, dlL 

2 It" " ( E i) E , s X i TT (B5) 


Then if equations (B2) and' (B3) are substituted into equa- 

^i ^ E i^R a a i ' 

tion (Bl) and -- ? ■■■ -r 1 is set equal to m , 

2CE i's,R 


dP t 

IT - Mi di + R 1+ l) - XiPi (B4) 
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APPENDIX C 


DEFINITION OF "e" TO A MATRIX POWER 


Inasmuch as [ A] is a matrix, e^ A ^ T must be defined. The 

([A]t) J 


symbol eW T is defined as 33 


3=0 


3 '- 


which may be evaluated 


more readily with the aid of the confluent form of Sylvester's 
Theorem (reference 3) . 


Let [u] be any square matrix of rank q with r distinct 
characteristic roots A 1 . . . X r , of multiplicity a x . . . o^., 

respectively. Then if P([u]) is a polynominal in [u] or a con- 
vergent infinite series in [u], the theorem states that 


r 

?([«]) = 

i=l 


1— 1 

X 
1 — 1 

V 

P(MTf (M - XkCE])" 

k*i 

\ 

(“» - l,! d ^ 1 

k. 

*** J 

V 

f 

x=\. 

> 


where [E] is the identity matrix of order q. 


(Cl) 


If all the roots are distinct then equation (Cl) reduces to 


P([u]) - S3 
i=l 


p(\) if (m -^m) 


k*i 


K) 


Example 1. - Calculate e when [u] = 


2 2 0 . Then the 

3 4 3 

characteristic roots are A 1 = 1, X 2 = 2, X 3 = 3, and 


10 0 
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»W = el 


(2[E] - [u])(3[E] - [u]) 

+ e 2 

(M - M)(3[E] - M) 

(2 - 1)(3 - 1) 


(1 - 2) (3 - 2) 

-J _ 


+ ©" 


([E] - [tt])(2[E] - M) 
(1 - 3)(2 - 3) 


from. which 


eW=a 


1 0 0 


0 0 0 


0 0 0 

-2 0 0 

+ e 2 

-2 -1 0 

+ e 3 

0 0 0 

2.5 0 0_ 


8 4 0_ 


_5.5 4 1_ 


Example 2. - Calculate eM where [u] '= 


2-2 3 

10 -4 5 

5-4 6 


Here A-^ = Ag = 1 and A^ = 2 so that 


5 [u] = 


_Z e A ([u] - X^Hg) 
dA A - 2 


A=1 


, e 2 ([u] -2^0 

(2- l) 2 



5 

4 

-8 


5 

2 

-5 


-4 

-4 

8 

= © 

15 

16 

-30 

+ © 

25. 

10 

-25 

+ e 2 

-15 

-15 

30 


10 

10 

-19 


15 

6 

-15_ 


.-10 

-10 

20. 



"lO 

6 

-13 


-4- 

-4 

‘ 8 

= © 

40 

26 

-55 

+ e 2 

-15 

-15 

30 


25 

16 

-34_ 


_-10 

-10 

20_ 
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